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Abstract It has been known since the 1960s that the layers of Na and K atoms, which occur between 80 and 
105 km in the Earth's atmosphere as a result of meteoric ablation, exhibit completely different seasonal behavior. 
In the extratropics Na varies annually, with a pronounced wintertime maximum and summertime minimum. 
However, K varies semiannually with a small summertime maximum and minima at the equinoxes. This contrasting 
behavior has never been satisfactorily explained. Here we use a combination of electronic structure and chemical 
kinetic rate theory to determine two key differences in the chemistries of K and Na. First, the neutralization of K + 
ions is only favored at low temperatures during summer. Second, cycling between K and its major neutral reservoir 
KHC0 3 is essentially temperature independent. A whole atmosphere model incorporating this new chemistry, 
together with a meteor input function, now correctly predicts the seasonal behavior of the K layer. 


1. Introduction 

The ablation of meteoroids in the upper atmosphere leads to the formation of distinct layers of metal atoms 
around 90 km. These layers, which are global in extent, provide a sensitive marker for the dynamics and 
chemistry of the atmosphere at the edge of space [Plane, 2003]. The presence of K atoms in the mesosphere 
and lower thermosphere (MLT) was first established from photometric observation of the twilight emission at 
766.5 and 769.9 nm [Sullivan and Hunten, 1 964]. Subsequent work using the resonance lidar technique [Megie 
eta!., 1978; Eska eta!., 1999] showed that the K column density varies semiannually, with maxima in summer and 
winter and minima at the equinoxes. The layer thus exhibits little overall seasonal variation at midlatitudes. In 
contrast, the Na layer has a pronounced annual variation, being around 3 times larger in winter than summer 
[Marsh eta!., 2013]. As a result, the Na/K abundance ratio at midlatitudes decreases from ~50 in winter to only ~10 
during summer. The explanation for this very surprising difference in behavior between two Group 1 (alkali) 
metals, whose chemistry should be very similar, is a problem of 50 years' standing. 

Initially, attention focused on whether the Na/K ratio indicated a seasonal change in the source of the metals in the 
MLT. The Na/K elemental ratio in meteoroids is around 1 5, whereas the seawater ratio is 47, suggesting that the 
vertical transport of sea-salt particles during storms in the polar night could be responsible for the wintertime 
enhancement of Na over K, which is greatest at high latitudes [Sullivan and Hunten, 1964; Megie et al., 1978], 
However, in addition to the challenge of transporting material upward within the winter polar vortex to the MLT 
where the mean wind is downward, Jegou et al. [1 980] showed that the mesospheric Na/Li ratio remains constant or 
may even decrease during winter, which appears to rule out a sea-salt source. An alternative explanation proposed 
that the neutral metal atom chemistry in the MLT is controlled by the dynamics of hydrated metal ion clusters [Jegou 
et al., 1 985], but the assumed temperature-dependent differences in the ion cluster chemistries required for this 
model have not subsequently been validated by experiment or theory. More recent modeling attempts [Eska eta!., 
1999; Delgado et al., 2006] have been unable to identify key differences in the chemistry of K and Na. 

K + ions have been observed by rocket-borne mass spectrometry in a permanent layer above 90 km [Kopp, 

1 997], These ions form continuously from K in the MLT via photoionization and charge transfer with ambient 
ions [Eska et al., 1999]: 


K + bv — > K + + e~ 

(1) 

K + NO+ -> K+ + NO 

(2) 

K + 0 2 + — > K + + 0 2 

(3) 
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Neutralization of K + occurs by forming a 
cluster ion with a ligand X, where X = N 2 , 0 2/ 
O (in order of decreasing concentration in 
the MLT). These clusters can then undergo 
ligand switching with another atmospheric 
component Y (C0 2 or H 2 0) to form a more 
stable cluster, before undergoing 
dissociative electron recombination: 

K+ + X(+ M) ^ K+ ■ X (4) 

K+-X + Y-+K+-Y + X (5) 
K + • (X or Y) + e~ — » K + X or V (6) 

where M is a third body (most likely N 2 or 0 2 ). 
Plane et at. [2006] carried out high-level 
quantum chemistry calculations on the K + -X 
and ■ X • K + • Y cluster ions and then used Rice- 
Ramsperger-Kassel-Markus (RRKM) theory to 
estimate the rate coefficients for reactions (4) 
and (5). The important result is that because 
K + is a relatively large singly charged ion (cf. 
Na + ), it forms very weakly bound clusters with 
the major species (N 2 , 0 2 , and O) in the MLT, 
with binding energies less than 20 Id mol -1 . 
The sequence of reactions (4)— (6) therefore 
only becomes significant at the very low 
temperatures characteristic of the 
summertime MLT. Figure 1 illustrates 
schematically the important pathways that 
couple K + and K in the MLT. 

Although the ion-molecule chemistry described above should produce an enhancement of K during summer, 
neutral chemistry would be expected, by analogy with Na [Plane, 2004], to remove K through formation of stable 
reservoir species. In the case of Na, the important reservoir is sodium bicarbonate (NaHC0 3 ), which is only 
converted back to atomic Na by photolysis or reaction with H. The reaction NaHC0 3 + H — > Na + H 2 C0 3 has an 
activation energy of about 10 kJ mor 1 [Cox et at., 2001], Thus, during summer the reaction becomes very slow, 
causing NaHC0 3 to build up and contributing to the observed minimum in the Na column density [Plane, 2004]. 

In this paper we use a combination of theoretical techniques to explore the analogous chemistry of 
potassium. The resulting rate coefficients for the neutral reactions, together with those from our previous 
study of K ion-molecule reactions [Plane et al., 2006], are then incorporated into a whole atmosphere 
chemistry-climate model. 

2. Theoretical Calculations on Neutral K Chemistry 

In order to calculate rate coefficients and photolysis rates, a variety of electronic structure calculations was 
performed using the Gaussian 09 suite of programs [Frisch eta!., 2009]. Initially, the hybrid density functional/ 
Flartree-Fock B3LYP method was employed with the 6-31 1 +G(2d,p) triple zeta basis set, which has both 
polarization and diffuse functions added to the atoms. Excited state calculations were performed using the 
time-dependent TD/B3LYP method [Bauernschmitt and Ahlrichs, 1 996]. Very accurate energies of the 
stationary points (i.e., local minima and transition states) on the relevant potential energy surfaces (PESs) 
were calculated at the complete basis set (CBS-QB3) level of theory [Montgomery et al., 1 999]. In the case of 
the KO radical, the CBS-QB3 method failed due to excessive mixing of the frozen core and valence orbitals in 
the coupled cluster step of the calculation. Accurate reaction enthalpies for reactions involving KO were 
therefore carried out at the B3LYP level using the augmented version of the relativistic small-core effective 
core potential ECP10MDF basis set of the Stuttgart/Cologne group [Lim et al., 2005]. 



Figure 1 . Schematic diagram of the important reaction pathways for 
K and K + produced by meteoric ablation in the MLT. Neutral and ionic 
K species are enclosed in green and blue boxes, respectively. Thin 
black lines depict chemical reactions; broken lines indicate photochemical 
processes. In this scheme, the main neutral reservoir is KFIC0 3 , which is 
permanently removed as a constituent of meteoric smoke. 


PLANE ETAL. 


©2014. The Authors. 


4754 




^AGU Geophysical Research Letters 


10.1 002/201 4GL060334 



Figure 2. (a) Potential energy surface for the recombination of KOH 
and C0 2 to form the reservoir species KHC0 3 . (b) Potential energy 
surface for the reaction between KHC0 3 and H, illustrating the three 
lowest-energy reaction channels. Color scheme: K (purple), O (red), 

H (white), and C (grey). 


As shown in Figure 1, K atoms are oxidized by 
0 3 to yield KO, which a second 0 3 can then 
oxidize to K0 2 [Plane, 2002]. K0 2 also forms 
directly from the recombination of K and 0 2 
[Plane et at., 1 990], Flowever, in the MLT these 
oxides exist in the presence of O, H, H 2 , Fl 2 0, 
and C0 2 . Calculations (see above) show that 
the following reactions are all exothermic: 

KO + O — > K + 0 2 

AH(0 K) = —243 kJ mol -1 (7) 

KO + H 2 0 -> KOH + OH 
AH(0 K) = -1 kJ mor 1 (8) 

KO + H 2 — > KOH + H 
AH(0 K) = -57 kJ mor 1 (9) 
K0 2 + O -> KO + 0 2 
AH(0 K) = -138 kJ mol -1 (10) 
K0 2 + H — > K 4- H0 2 
AH(0 K) = -85 kJ mor 1 (11) 
KOH + H K + H 2 0 
AH(0 K) = -171 kJ mor 1 (12) 

Reaction (9) has a second channel forming 
K + H 2 0 which is exothermic by 
228 kJ mor 1 . However, in the case of the 
analogous NaO + H 2 reaction, formation of 
Na is the minor channel [Ager and Howard, 
1987] and we assume the same for K. 
Reaction (11) could also form KOH + O, which 
is exothermic by 135 kJ mor 1 but involves 
breaking the strong O — O bond in the 
superoxide and is therefore probably slow. 


The potential energy surfaces of reactions (7)— (1 2) do not contain significant energy barriers, and so by analogy 
with the corresponding reactions of NaO and Na0 2 [Ager and Howard, 1 987; Helmer and Plane, 1 993; Griffin 
et at., 2001], these reactions are likely to be fast with little temperature dependence. We therefore assign 
them a typical collision theory rate coefficient of 2 x 10 -10 exp(— 120/7) cm 3 molecule -1 s 


KOH can also recombine with C0 2 to form the bicarbonate: 

KOH + C0 2 (+ M) — > KHC0 3 AH(0 K) = - 1 42 kJ mol -1 (1 3) 

Figure 2a illustrates the potential curve for reaction (13), together with the geometries of the stationary 
points (atom coordinates are listed in Table SI in the supporting information). The KOH and C0 2 initially 
recombine to form the K(OH)C0 2 adduct. This adduct then rearranges to the more stable KHC0 3 either 
via migration of the H or the K atom, the latter having a much lower energy barrier with respect to the 
energy of the reactants. The rate coefficient for reaction (1 3) was calculated using a multiwell energy-grained 
master equation solved with the program Master Equation Solver for Multienergy Well Reactions 
[Robertson et at., 2008], where the microcanonical rate coefficients linking the reactants, intermediates, and 
products were calculated by RRKM theory using the molecular parameters in Table SI (supporting 
information). The internal energies of the intermediates on the PES were divided into a contiguous set of 
grains (width 100 cm -1 ), each containing a bundle of rovibrational states. For dissociation to products or 
reactants, microcanonical rate coefficients were determined using inverse Laplace transformation to link 
them directly to the dipole-induced dipole capture rate coefficient [Georgievskii and Klippenstein, 2005], 
kcapture(KOH + C0 2 ) = 5.8 x 10 -10 (7/ 300 K) 1/6 cm 3 molecule -1 s -1 . The probability of collisional transfer 
between grains was estimated using the exponential down model (average energy for downward transitions 
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was set to < A£ > down = 300 cm \ assumed temperature independent). At the low-pressure limit appropriate 
for the MLT region, /c 13 (1 40-220 K) = 7.1 x 10 -28 (200 K/T ) 4 ' 21 cm 6 molecule -2 s -1 . 

KHCO3 could potentially be destroyed by reaction with H or O: 

KHCO 3 + H — ► K + H 2 C0 3 AH(0 K) = 1 Id mol -1 (14a) 

■^KC0 3 + H 2 AH(0 K) = -3 Id moP 1 (14b) 

^K+H 2 0 + C0 2 AH(0 K) = -29 kJ mol -1 (14c) 

KHCO 3 + O-+KCO 3 + OH AH(0 K) = 8 kJ mol -1 (15) 

These reactions are either exothermic or only slightly endothermic and so were investigated in more detail. 
Figure 2b shows the potential energy curves and transition state geometries for the three channels of 
reactions (1 4a)— (1 4c). Channels (1 4a)— (1 4c) have barriers of 32, 61, and 71 kJ mol -1 , respectively. Channels 
(14b) and (14c) can therefore be ruled out at MLT temperatures. The rate coefficient for reaction (14a) was 
estimated to be /c 14a (1 40-220 K) = 4.5 x10 -11 exp(— 3590/T) cm 3 molecule -1 s -1 by using transition state 
theory with a correction for tunneling of FI through the reaction barrier, assumed to be a symmetrical 
Eckart barrier as in our previous study of the analogous reaction of NaHC0 3 [Coxetal., 2001] (see Table SI in 
the supporting information for the relevant molecular parameters). At the various levels of theory used 
in this study, the energy barrier for KFICO3 + FI is consistently between 7.1 and 9.3 kJ mol -1 higher 
than that for NaHC0 3 + FI. Thus, even allowing for some uncertainty in the calculated absolute value of 
the barrier for reaction (14a), the reaction should be at least 2 orders of magnitude slower at 200 K 
than its Na analogue. Reaction (15) has a significant barrier of 38 Id mor 1 and will therefore be too 
slow in the MLT for further consideration. 

The photolysis cross sections for KHC0 3 , K0 2 , and KOH as a function of wavelength were calculated at the TD/ 
B3LYP/6-31 1 +g(2d,p) level of theory [Frisch et al., 2009] and then placed on an absolute scale using our 
previously measured cross sections for the Na analogues [Self and Plane, 2002]. The photolysis rates of these 
three potential molecular reservoirs were estimated to be 1.2 x 10 -4 , 2.2 x 1 0 3 , and 0.027 s -1 , respectively. 
Figure SI (supporting information) illustrates the cross-section spectra for KFICO 3 and NaHC0 3 , which have 
thresholds around 280 nm. 

Finally, the dimerization rate coefficient of KFICO 3 was estimated from the dipole-dipole capture rate 
[Georgievskii and Klippenstein, 2005] to be 1.0 x 10 -9 (7/ 200 K) - 029 cm 3 molecule -1 s -1 , using a calculated 
value for the dipole moment of KFICO 3 of 8.2 Debye (at the B3LYP/6-31 1 + g(2d,p) level). Plowever, KPICO 3 can 
also polymerize with other meteoric constituent molecules (e.g., NaFHC0 3 , FeOFI, and Mg(OFI) 2 ), and the 
dimerization rate coefficient needs to be increased to account for this. The meteoritic elemental ratios of Na, 
Fe, Mg, and Si relative to K are 1 5.5, 248, 279, and 264, respectively [Asplund et al., 2009]. Thus, the 
concentration of compounds containing these other elements could potentially be -800 times higher than 
the concentration of KFIC0 3 . Flowever, recent modeling studies [Vondrak et al., 2008; Feng et al., 2013; 
Langowski et al., 201 4] indicate that Fe, Mg, and Si ablate with much lower efficiencies than K (and Na), so the 
polymerization rate coefficient should be estimated by increasing the dimerization rate coefficient by 
significantly less than 800. A factor of 270 (i.e., a phenomenological polymerization rate coefficient of 
2.7 x 10 - 7 cm 3 molecule -1 s -1 ) was found to reproduce the observed underside of the K layer satisfactorily, 
as shown below. 

3. Atmospheric Modeling 

Table 1 contains a complete list of rate coefficients for the neutral and ion-molecule reactions depicted in 
Figure 1. In order to model the K layer, we adopt the same approach that was used recently to model the 
global Na, Fe, and Mg layers [Feng et al., 2013; Marsh et al., 2013; Langowski et al., 2014]. Three components 
are combined to construct the K model (termed here WACCM-K): (i) the Whole Atmosphere Climate 
Community Model (WACCM, version 4), a "high-top" coupled chemistry-climate model with an upper 
boundary at 6.0 x 10 -6 hPa (-140 km); (ii) a description of the neutral and ion-molecule chemistry of K 
(Figure 1 and Table 1); and (iii) the meteoroid input function (MIF), which specifies the temporal and 
geographical input of meteoric metallic atoms into the MLT [Marsh et al., 2013]. 
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Table 1 . Neutral and Ionic Gas-Phase Reactions in the Potassium Model 

Number 3 Reaction Rate Coefficient 13 Source 


Neutral Chemistry 


R1 

k+o 3 ->ko + o 2 

1.15X 10 -9 exp(— 120/T) 

1 

R2 (7) 

ko + o^k+o 2 

2 x 1 0 - 1 0 exp(— 1 20/7) 

2 

R3 

k+o 2 (+ M)->K0 2 

1.3 x10 -29 (77200) - 1 23 

3 

R4 

ko + o 3 -»ko 2 + o 2 

6.9x 10 -1 ° (exp(— 385/7) 

1 

R6 (10) 

ko 2 + o->ko + o 2 

2x 10 - 10 exp(- 120/7) 

2 

R7 (8) 

ko + h 2 o~>koh + oh 

2x 10 - 10 exp(- 120/7) 

2 

R8 (9) 

ko + h 2 ->koh + h 

2x 10 - 10 exp(- 120/7) 

2 

R9 (11) 

ko 2 + h^k+ho 2 

2 x 1 0 - 1 0 exp(— 1 20/7) 

2 

R10 (12) 

koh + h^k+h 2 o 

2 x 1 0 - 1 0 exp(— 1 20/7) 

2 

R11 (13) 

koh+co 2 (+ M)^KHC0 3 

7.1 xIO -28 (T/200) 4 2 

4 

R12 (14a, 14b, 14c) 

khco 3 + h^k+h 2 co 3 

4.5 x10 -11 exp(— 3590/7) 

5 

R13 

KHC0 3 + KHC0 3 (+ M) -> dimer 

Ion-molecule Chemistry 

270 x [1 X 10 9 (T/200) °- 29 ] 

6 

R20 (2) 

K + NO + — > K + + NO 

9.4 X 10 — 10 

7 

R21 (3) 

k+o 2 + ->k + + o 2 

3.2 xIO -9 

7 

R22 (4) 

K + + N 2 (+ M)->K + N 2 

2.3 xIO -30 (77200) - 239 

7 

R-22 

K + N 2 (+M)-*K + + n 2 

2.8 xl0 -8 exp(- 1680/7) 

7 

R23 (4) 

K + + 0 2 (+ M)->K + 0 2 

1 .2 x 1 0 — 30 ( 77200) 2-1 2 

7 

R-23 

K + 0 2 + (+ M)-»K + + 0 2 

1 .5 x 1 0 -9 exp(— 820/7) 

7 

R24 (4) 

K + + O (+ M) -> K + ■ 0 + 

8.8 xIO -32 (7/200 K) -1 - 28 

7 

R-24 

K + ■ 0 + (+ M) — > K + + O 

2.6 x10 -1 °exp(-1 800/7) 

7 

R25 (4) 

K + + C0 2 (+ M)-»K + C0 2 

1.3 xIO -29 (7/200) - 243 

7 

R26 (4) 

K + + H 2 0 (+ M)->K + H 2 0 

3.0 xIO -29 (7/200) - 222 

7 

R27 (5) 

k + n 2 + o->k + o + n 2 

2.9 X 1 0 10 (7/200 K) -0 - 17 

7 

R-27 (5) 

k + o + n 2 ->k + n 2 + o 

2.5 xIO -11 (7/200 K) -0 - 55 

7 

R28 (5) 

k + n 2 + co 2 ->k + co 2 + n 2 

4.8 X 1 0 10 (7/200 K) -0 - 88 

7 

R-28 (5) 

k + co 2 + n 2 -*k + n 2 + co 2 

2.8 x 10 -1 ° exp(— 2220/7) 

7 

R29 (5) 

k + o 2 + o^k + o + o 2 

2.8 xIO -10 (7/200 K) -0 - 42 

7 

R-29 (5) 

k + o + o 2 ^k + o 2 + o 

5.0x 10 -1 ° exp(— 752/7) 

7 

R30 (5) 

k + o+co 2 ^k + co 2 + o 

7.1 xIO -10 (7/200 K) -0 - 21 

7 

R-30 (5) 

k + co 2 + o->k + o + co 2 

1 .4 x 1 0 -9 exp(— 2200/7) 

7 

R31 (5) 

k + o + h 2 o->k + h 2 o+o 

7.1 xIO -10 (7/200 K) -1 - 90 

7 

R32 (5) 

k + o 2 + n 2 ->k + n 2 + o 2 

1.6x 10 10 (7/200 K) -0 - 90 

7 

R-32 (5) 

K + ■ N 2 + 0 2 — > K + ■ 0 2 + N 2 

I.Ox 10 -9 exp (-873/7) 

7 

R33 (5) 

k + o 2 + co 2 ^k + co 2 + o 2 

1.5 x 10 - 10 (7/200 K) -1 - 94 

7 

R34 (5) 

k + o 2 + h 2 o->k + h 2 o+o 2 

1.8 xIO -9 (7/200 K) -0 " 

7 

R35 (5) 

k + n 2 + h 2 o->k + h 2 o + n 2 

2.4 xIO -9 (7/200 K) -045 

7 

R36 (5) 

K + ■ C0 2 + H 2 0 — > K + ■ H 2 0 + C0 2 

1.4X10 -9 (7/200 K) -1 - 26 

7 

R37 (6) 

KX + + e - — > K + X (X = O, 0 2: N 2 , C0 2 , H 2 0) 
Photochemical Reactions 

1 xIO -6 (7/200) 1 72 

8 

R38 

K0 2 + /»'—> K + 0 2 

2.2x 10 3 

9 

R39 

KOPI + /)v-»K + OH 

2.7 X 10 2 

9 

R40 

KHC0 3 + hv-»K+HC0 3 

1.2X10 -4 

9 

R41 (1) 

K + hv^K + + e - 

4x10 -5 

10 


a The number in parenthesis refers to the number used when the reaction is listed in the text. 
b Rate coefficient units: unimolecular, s -1 ; bimolecular, cm 3 molecule -1 s -1 ; and termolecular, cm 6 molecule -2 s -1 . 
1 . Measured [Plane, 2002], 2. Estimated rate coefficient from electronic structure calculations on the relevant transition 
states (see text). 3. Measured [Plane et al. : 1 990], 4. Electronic structure calculations and RRKM theory (see text). 5. 
Electronic structure calculations and TST (see main text). 6. This reaction represents loss of KHC0 3 to meteoric smoke 
particles (see text). 7. Theory [Plane et at., 2006], 8. Estimate based on typical dissociative electron recombination 
reaction rates [Florescu-Mitchell and Mitchell, 2006], 9. Electronic structure calculations using time-dependent 
density function theory (see text). 10. Estimate [Swider, 1970], 


WACCM-K was run from 2004 to 2012 (after a 9 year spin-up), with dynamical fields nudged in the 
troposphere and stratosphere to the Goddard Earth Observing System-5 (GEOS-5) meteorological data set, as 
described by Feng et al. [2013]. A nudging coefficient value (0.01) was used when assimilating the GEOS-5 
analysis into WACOM so that 1 % of the meteorological conditions were combined with WACCM fields below 
60 km at every model dynamics time step. Above 60 km the model was free running. The MIF is calculated 
from current knowledge of the astronomical characteristics of the sporadic meteor complex to estimate the 
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The WACCM-K output was compared 
with lidar measurements of the K layer 
made over 3 years (2004-2006) at 
Kuhlungsborn, Germany (54°N, 12°E) 
(see Eska et at. [1 999] for a description 
of the technique). The K layer at this 
station is typical of other midlatitude 
locations [Sullivan and Hunten, 1964; 
Megie et at., 1978]. To minimize the 
impact of solar cycle effects on the K 
layer, we focus here on the seasonal 
variation of the K layer averaged over 
the 3 years. Figure 3a shows the 
measured seasonal variation of the K 
density as a function of height, and 
Figure 3b shows the modeled 
variation. The first noteworthy point is 
that the semiannual variation of the K 
layer, with a midsummer maximum 
and minima at the equinoxes, is very 
well captured by the model. A second 
point is that the absolute K density 
agrees well with the observed density, 
using a MIF for K (average 120 atom 
cm _2 s _1 ) which is in the chondritic 
ratio (K/Na = 6.5% [Asplund et al., 
2009]) to the Na MIF used in our recent 
study of the Na layer [Marsh et al., 

201 3]. This is in accord with a chemical 
ablation model [Vondrak et al., 2008] 

which predicts that Na and K should ablate with similarly high efficiencies compared to other meteoric 
metals. A third point is that the peak height and width of the layer are satisfactorily captured by the model 
(see Figure S2 in the supporting information which shows MLT profiles of all the K-containing species in 
WACCM-K). Figure 3c shows that the modeled semiannual seasonal variation of the K column abundance is 
also in good agreement with the lidar measurements. The contrasting seasonal variation of the Na layer at 
the same location [Marsh et al., 2013] is included for comparison. 



8 E 

o 

7 I 

flj 

6 


Figure 3. Annual variation of the K layer at Kuhlungsborn, Germany (54°N, 
12°E), averaged from 2004 to 2006. (a) K density (atom crrT 3 ) as a function of 
altitude and month, measured by lidar. (b) K density as a function of altitude and 
month from WACCM-K. (c) The observed and modeled K column abundance, 
compared with the Na column abundance (since Na lidar measurements are 
not available from Kuhlungsborn for the same period, validated model output 
[Marsh et al., 2013] for that location averaged over 2004-2006 is shown). 


global meteoric mass flux deposited at 
each altitude within the MLT at all 
latitudes for each day of the year 
[Fentzke and Janches, 2008]. A 
chemical ablation model [Vondrak 
et al., 2008] is then used to compute 
the injection rates of individual 
metallic elements as a function of 
altitude. See Feng et al. [2013] for 
further details. 


4. Discussion and Conclusions 

Why does K exhibit a summertime peak and equinoctial minima, in contrast to Na? There are two factors 
involved. First, the K + ion is more difficult to neutralize because it can only form weakly bound cluster ions at 
very low temperatures. Thus, the low temperatures of the summertime MLT shift the balance from K + to K, 
accounting for the summertime maximum. The second factor concerns the neutral chemistry. None of the 
reactions involving KO, K0 2 , and KOFI has a significant temperature dependence. In contrast, the activation 
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energy for the reaction KHC0 3 + H is so large that this reaction never plays a significant role, even during the 
warmest period in winter. The only route from KHC0 3 back to K is via photolysis, and this is essentially 
temperature independent. For Na, the reaction NaFIC0 3 + FI has a smaller activation energy [Coxetal., 2001], 
and so during winter this reaction becomes fast enough to recycle the reservoir to atomic Na, leading to the 
pronounced wintertime maximum. Because there is less winter-summer variation in the K layer, the 
semiannual variations become pronounced. As shown in Figure S3 (supporting information), atomic O 
exhibits minima at the equinoxes, whereas Fl 2 and C0 2 exhibit maxima. Inspection of Figure 1 shows that 
simultaneously decreasing O and increasing Fl 2 and C0 2 will drive K to KFIC0 3 , and this explains the minima 
in the K layer seen during the equinoxes. 

It should be pointed out that the MIF used here would correspond to a global input rate of 4.6 1 cT 1 of 
ablated cosmic dust [ Marsh eta!., 201 3]. This is at the lower end of current estimates of this quantity, which vary 
from ~3 to 300 1 cT 1 [Plane, 2012]. Although this estimate is calculated from an astronomical model constrained 
to measurements of meteor head echoes using a number of high-performance large aperture radars 
[Fentzke and Janches, 2008], it is possible that a population of undetected meteors could make a significant 
contribution to the MIF. To test this possibility, we ran a sensitivity experiment where the MIF was increased by a 
factor of 1 0 in WACCM-K. The K column density increased by a factor of 8.6 (there is a small nonlinearity 
because of the second-order polymerization reaction R13 in Table 1). The coefficient of R13 would therefore 
need to be increased if the MIF was larger than 4.6 1 d _1 . Flowever, the layer height, the layer shape, and the 
semiannual seasonal behavior of the K layer — which is the focus of this study — were unchanged. 

In conclusion, the surprisingly different behavior of the K layer now appears to be explained. Nevertheless, it 
would be desirable if the theoretical rate calculations presented here for reactions (13) and (14a)-(14c), and 
the assumed rate coefficients for reactions (8)— (1 2), could be tested experimentally in the future. 
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